Variational methods, multiprecision and nonrelativistic energies 
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It is known that the variational methods are the most powerful tool for studying the Coulomb 
three-body bound state problem. However, they often suffer from loss of stability when the number 
of basis functions increases. This problem can be cured by applying the multiprecision package 
designed by D.H. Bailey. We consider the variational basis functions of the type exp( — a n ri —/3 n r2 — 
7n?*i2) with complex exponents. The method yields the best available energies for the ground states 
of the helium atom and the positive hydrogen ion as well as many other known atomic and molecular 
systems. 



1. The development of the variational method for the 
Coulomb bound state problem can be traced using as an 
example the ground state of the helium atom. In early 
days when computers were big and very expensive the 
search proceeded mainly in the direction of making ex- 
pansion of the variational wave function as compact as 
possible (in a sense of number of variational parameters 
and/or basis sets). At first, the explicitly correlated basis 
were introduced [00] now called as the Hyllcraas basis 



s = n + r 2 , u = ru, t = + r 2 , 

then it became clear that at least for the ground state 
of the helium atom it is essential to incorporate into 
the wave function such peculiarity as the logarithmic be- 



haviour of the type RhxR at R 



„2U 



0, 



first analytically derived independently by Bartlett and 
Fock j|. In 1966, Frankowski and Pekeris (see Table ||) 
introduced the compact representation W of the form 



^(ri,r 2 )=e K 



(,s 2 + t 2 y/ 2 (in s y 



and later, in 1984, Freund and co-workers || reported 
even more compact expansion of the same form. Inclu- 
sion of the logarithmic term into the variational wave 
function brought substantial improvement of nonrela- 
tivistic energies for the two electron atoms. In 1994, 
Thakkar and Koga || have found a compact expansion 
without logarithms which uses powers that are not inte- 
gers nor even half integers. As far as we know none of 
these compact expansions has been used for analytical 
evaluation of matrix elements of the Breit interaction. 

With advance of computer power basis sets became 
simplified that allowed for calculation of numerous ma- 
trix elements required for relativistic and QED correc- 
tions. The efforts were concentrated on a choice of 
a strategy that defines a sequence of basis functions 
genereated. In the double basis set method with gen- 
eralyzed Hylleraas basis functions 



V(ri,r 2 ) = J2^Hri2e- ari - 0r2 

+ / , Ctjfc' 1' 2 r 12 e 

were used. This double basis set technique along with 
full optimization of nonlinear parameters at each basis 
set yield substantial progress in accuracy. However, the 
main factor that hinder further advance become the nu- 
merical instability due to almost linear dependence of the 
basis set at large N . 

The work of Goldman || is a bit apart of the main 
path. It recovers the idea of Pekeris Q| to use uncoupled 
coordinates and orthogonal Laguerre and Jacoby polyno- 
mials as basis functions. 

The method expounded in our work is a continuation 
of efforts by Drake and Yan to utilize as much simple 
basis functions (geminals) as possible. 

2. Expansion we want to consider here is very similar 
to the generalized Hyllcraas basis set, but instead of us- 
ing the polynomials over Hyllcraas variables we generate 
nonlinear parameters in the exponents in a quasi-random 
manner, 



-ari— f3r 2 — 77-12 



-air 1 —Pir2—~/iri2 



(1) 



This method has been successfully used in calculations 
P,^0[ previously. Obviously, the matrix elements can be 
evaluated in the same way as for the generalized Hyller- 
aas basis set (Q). Moreover, if one replaces real exponents 
by complex exponents the integrals will remain exactly 
the same as for the real case. In its strategy the method 
is very close to the SVM method by Varga, Suzuki [pd]| , 
where gaussians are exploited instead. 

In a formal way, a variational wave function is ex- 
panded in a form 

00 

tpo = \Ui _Re[exp {-ctin ~ f5 t r 2 - 7^12)] 
1=1 

+W l Im[exp{~a i r 1 -(3 l r 2 - ll r 12 )]}y l L i ^(r 1 ,r 2 )- 

Here cti, f3i and 7^ are complex parameters generated in 
a quasi-random manner 1 1 3 14J]: 



1 



Oii 



i(i + l)^/pZ 



[(A 2 - Ax) + A x ]+ 
(A'.-A'J + AW 



[x\ designates the fractional part of x, p a and q a are 
some prime numbers, [Ai,^] and [yl^,^] are real vari- 
ational intervals which need to be optimized. Parameters 
f3i and 7, are obtained in a similar way. 

An important feature of the method is that it demon- 
strates a very fast convergence. The general rule which 
can be inferred experimentally from the use of the 
method is that increasing of the basis by about 200 func- 
tions yields about one additional digit in the variational 
energy. The minor deficiency is that the basis quickly 
degenerates when TV increases. Already for moderate 
TV ~ 250 — 400 a quadruple precision is required. 

Multiprecision package of Fortran routines MPFUN 
has been designed by David H. Bailey for computa- 
tions with floating point numbers of an arbitrary length. 
Usually it is necessary to make significant changes into 
Fortran source code in case if Fortran-77 language is used. 
Fortunately, the author of MPFUN package has devel- 
oped a translator program that facilitate converting the 
programs to multiprecision drastically. In general, two 
directives incorporated as comments in a source code are 
required per one routine. For example a source code for 
the considered variational method has been transformed 
to multiprecision version within two hours of manual 
work. Eventually a code we've gotten has been tested 
on a personal computer with the Celeron 500 MHz pro- 
cessor. For one run with the basis of TV — 1400 functions 
and 40 decimal digits it requires about 3 hours. 

For users of Fortran-90 no preprocessor is needed due 
to new advanced features of Fortran-90, such as derived 
data types and operator extensions. 



TV 


E (a.u.) 


1400 


-2.90372437703411959629 


1600 


-2.903724377034119597843 


1800 


-2.9037243770341195981964 


2000 


-2.9037243770341195982713 


2200 


-2.9037243770341195982955 


extrapolation 


-2.903724377034119598306(10) 



and TV = 1600. Quadruple precision was not sufficient 
at these TV and we used the multiprecision version of the 
program with 40 significant decimal digits. Further cal- 
culations with TV = 1800 — 2200 were performed with 
48 significant digits and only partial optimization of the 
parameters of the last layer (corresponding to the region 
where the logarithmic behaviour is the most essential) 
was done. Some optimization of a distribution of rij be- 
tween the layers (TV = m + ri2 + U3 + n^) was carried out 
as well. 

As can be seen from the Tabic [rj the present result 
extends the accuracy of the nonrelativistic ground state 
energy for the helium atom by as much as 3 decimal dig- 
its. 



N 



E (a.u. 



Frankowski and 

Pekeris [§ 
Freund, Huxtable 

and Morgan III 
Thakkar and Koga 
Drake and Yan 
Goldman [| 
This work 



246 -2.9037243770326 

230 -2.9037243770340 

308 -2.9037243770341144 

1262 -2.90372437703411948 

8066 -2.903724377034119594 

2200 -2.903724377034119598296 



TABLE II. Comparison of the ground state energy of the 
helium atom obtained in this work with other theoretical cal- 
culations. 

Second case is the hydrogen molecular ion ground state 
that represent an other limit of mass distribution of con- 
stituents with one light and two heavy particles. For 
this case it is especially essential that we introduce com- 
plex exponents, because it is the most natural way to 
suit the oscillatory behaviour of the vibrational motion 



in the wave function. In this case (see Table [II) again 
40 decimal digits have been used for TV — 1400 — 1800 
and 48 decimal digits for large TV to provide the numer- 
ical stability of the calculations. Table IV demonstrates 
progress in obtaining variational nonrelativistic energy 
for this state. The accuracy is extended by as much as 4 
additional digits. 



TABLE I. Variational energy (in a.u.) of the helium 
ground state as a function of TV, the number of basis func- 
tions. 

In our calculations for the helium ground state four ba- 
sis sets with independently optimized nonlinear param- 
eters were used. These sets were built up like a pine 
tree. The first layer was tuned to approximate the gen- 
eral behaviour of the solution at intermediate and large 
r\ and T2 ■ The second layer was chosen to be flexible in a 
smaller region of r\ and r-i and so forth. A detailed opti- 
mization was performed for the sets with total TV = 1400 



TV 


E (a.u.) 


1400 


-0.597139063123404975 


1600 


-0.597139063123405047 


1800 


-0.5971390631234050655 


2000 


-0.5971390631234050710 


2200 


-0.5971390631234050740 


extrapolation 


-0.597139063123405076(2) 



TABLE III. Variational energy (in a.u.) of the positive 
hydrogen ion ground state as a function of TV, the number of 
basis functions. 



2 



N 



E (a.u.; 



Gremaud, Delande 
and Billy |l| 
Rebane and Filinsky 
Moss |l7| 
This work 



31746 -0.597139063123 



2200 



-0.59713906312340 

-0.5971390631234 

-0.597139063123405074 



TABLE IV. Comparison of the ground state energy of the 
positive hydrogen molecular ion obtained in this work with 
other theoretical calculations. m v — 1836.152701m e . 

In Table |y| the other examples are summarized. A 
negative positronium ion demonstrates a limit of three 
particles of equal masses. The second and third cases are 
applications of the method to the states with nonzero 
angular momentum. The last example in this Table is 
of special proud. That is the last vibrational state in a 
series of 5*-states of the hydrogen molecular cation, and 
that is the first variational confirmation of the existence 
of this state (the binding energy corresponding to the 
cited value is 0.74421(2) cm -1 ). The accuracy of the 
artificial channels scattering method (2^] is presumably 
better, however, wave functions are not forthcoming with 
this method that makes difficult calculation of physical 
properties of the state other than energy. 



system 



E 



He(2 3 P) 



4 He+p(L = 35,u = 0) This work 



H+ (L = 0,u = 19) 



This work -0.2620050702329801077(3) 
H -0.262005070232976 

This work -2.13316419077928310(2) 
[[| -2.13316419077927(1) 

-2.98402095449725(1) 
-2.98402094 

This work -0.4997312306 
[fn) -0.49973123063 



TABLE V. Other examples of three-body calculations. (L 
is the total angular momentum, v is the vibrational quantum 
number.) 

3. One may say that this high accuracy is redundant 
and has no physical meaning. But obviously, it shows the 
power of modern computers and theirs ability to solve the 
quantum three-body problem to any required accuracy. 
On the other hand, uncertainty in the variational wave 
function approximately as much as the square root of 
the uncertainty in the variational energy and is about 
10~ 9 — 1CP 10 . This accuracy does not look redundant. 
These results prove that the nonrelativistic bound state 
three-body problem is now satisfactorily solved and the 
main efforts should be addressed to relativistic and QED 
effects. 

The other advantage of the method is the simplicity of 
the basis functions that allows for evaluate analytically 



relativistic matrix elements of the Breit Hamiltonian. It 
is possible as well to evaluate analytically the vacuum 
polarization term (Uehling potential) |||] and to build 
up an effective numerical scheme for the one-loop self- 
energy corrections |23| . These features make the consid- 
ered variational method to be highly powerful universal 
tool for studying the three-body problem. 

This work has been partially supported by INTAS 
Grant No. 97-11032, which is gratefully acknowledged. 
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